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ABSTRACT 

We present the analysis of a Chandra High-Energy Transmission Grating (HETG) 
observation of the local Seyfert galaxy NGG 1365. The source, well known for its 
dramatic X-ray spectral variability, was caught in a reflection-dominated, Compton- 
thick state. The high spatial resolution afforded by Chandra allowed us to isolate 
the soft X-ray emission from the active nucleus, neglecting most of the contribution 
from the kpc-scale starburst ring. The HETG spectra thus revealed a wealth of He- 
and H-like lines from photoionized gas, whereas in larger aperture observations these 
are almost exclusively produced through collisional ionization in the circumnuclear 
environment. Once the residual thermal component is accounted for, the emission¬ 
line properties of the photoionized region close to the hard X-ray continuum source 
indicate that NGC 1365 has some similarities to the local population of obscured active 
galaxies. In spite of the limited overall data quality, several soft X-ray lines seem to 
have fairly broad profiles (~800-1300 km s“^ full-width at half maximum), and a range 
of outflow velocities (up to ^1600 km s“^, but possibly reaching a few thousands km 
s“^) appears to be involved. At higher energies, the Ka fluorescence line from neutral 
iron is resolved with > 99 per cent confidence, and its width of ~3000 km s“^ points 
to an origin from the same broad-line region clouds responsible for eclipsing the X-ray 
source and likely shielding the narrow-line region. 

Key -words: galaxies: active - X-rays: galaxies - galaxies: individual: NGC 1365 - 
line: identification 


1 INTRODUCTION 

Irrespective of their different optical properties and classi¬ 
fication, the energy source of active galactic nuclei (AGN) 
is strongly suggested to be intrinsically the same (e.g. An- 
tonucci & Miller 1985). When broad lines are lacking, a di¬ 
rect view of central engine would be prevented by a dusty 
absorber, covering a substantial, yet not complete fraction of 
the solid angle, and thus commonly envisaged as a toroidal 
structure (Krolik & Begelman 1988). The orientation-based 
idea that underlies AGN unification models (Antonucci 
1993, and references therein) is still broadly accepted. While 
it was argued early on that the torus might well consist of 
individual clouds, the multiple pieces of observational evi¬ 
dence accumulated in the last decade imply that the cir¬ 
cumnuclear environment in AGN is characterized by a large 
degree of dumpiness and inhomogeneity on various physical 
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scales, from tens of pc down to a few hundreds of gravita¬ 
tional radii (rg = GMbh/c^) from the supermassive black 
hole (see Bianchi, Maiolino & Risaliti 2012; Netzer 2015 for 
recent overviews). 

The nearby Seyfert galaxy NGC 1365 {z ~ 0.0055) is an 
ideal object to probe the complexity of nuclear obscuration 
in AGN. Optically classified as a type 1.8 (Maiolino & Rieke 
1995), NGC 1365 has been extensively studied in the X- 
rays for its exceptional variability, which makes it the most 
striking example of a changing-look AGN (Matt, Guainazzi 
& Maiolino 2003). Apparently due to a favourable line of 
sight, the inner accretion disc/X-ray corona system is seen 
through an erratic absorbing medium, whose marked col¬ 
umn density gradients are betrayed by the random alterna¬ 
tion of Compton-thin (Ah < 10^^ cm“^) and Compton-thick 
(Ah > 10^^ cm“^) states. The switch from transmission- to 
reflection-dominated spectra (and back) can take less than 
two days, involving the transit of a single cloud (Risaliti et 
al. 2007). This kind of eclipses not only sets a tight empir- 
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ical limit to the size of the X-ray source (no larger than a 
few tens of Vg), but also provides an implicit measure of the 
distance and density of the intervening blobs (both indica¬ 
tive of the gas in the Broad Line Region; BLR), and reveals 
unique details about their shape (e.g. Maiolino et al. 2010). 
Besides this, in the Compton-thin regime NGC 1365 regu¬ 
larly exhibits four strong absorption lines in the ~6.7-8.3 
keV band. Their spacing prompts an identification with the 
Ka and K/3 pairs from Fexxv-xxvi, blueshifted by reason 
of an outflow velocity of Uout ~ 1000-5000 km s“^, vari¬ 
able over timescales of weeks/months (Risaliti et al. 2005; 
Nardini et ah, in preparation). These lines could just trace 
the high-ionization phase of a much more pervasive wind, 
possibly including all the manifest absorption components 
(Connolly, McHardy & Dwelly 2014; Braito et al. 2014). 

Adding to the overall merit of NGC 1365 is one of the 
largest infrared luminosities (Lir ~ 10^^ ^q) among nearby, 
non-interacting galaxies. This mostly arises from the great 
number of star-forming clusters in the central ~2-3 kpc, 
disposed in an elongated ring around the nucleus (Alonso- 
Herrero et al. 2012, and references therein). The powerful 
starburst is also the main source of the observed soft X- 
ray emission (E < 3 keV), which is found to be constant 
in flux and shape. Indeed, the ~500-ks high-resolution spec¬ 
trum obtained by combining all the 2004 and 2007 data from 
the XMM-Newton Reflection Grating Spectrometer (RGS) 
is rife with collisionally ionized lines from the hot diffuse 
gas, at odds with a typical obscured Seyfert (Guainazzi et 
al. 2009). In this paper we report on two Chandra observa¬ 
tions of NGC 1365 taken with the High-Energy Transmission 
Grating (HETG; Canizares et al. 2005) over a span of four 
days, during which the AGN remained in a heavily absorbed, 
Gompton-thick state. 


2 OBSERVATIONS AND DATA REDUCTION 

NGG 1365 was observed twice by Chandra in 2012, on April 
9 and three days later (ObsIDs 13920/1), for a total expo¬ 
sure of ~200 ks. After a series of six 15-ks snapshots in 2006 
(Risaliti et al. 2007) and a previous one of the same length in 
2002, these were the hrst observations of NGC 1365 at high 
spectral resolution with Chandra. The HETG consists of two 
sets of gratings, optimized for medium (MEG; 0.5-7 keV) 
and high (HEG; 0.8-10 keV) energies, respectively, and was 
used in combination with the Advanced CGD Imaging Spec¬ 
trometer (AGIS-S) array. The data were reprocessed with 
the CIAO version 4.4 software package and the v4.4.9 Cali¬ 
bration Database (CALDB). For each observation, the MEG 
and HEG spectra were extracted from hrst-order events in 
both (±1) diffraction arms. Redistribution matrices were 
created with the CIAO tool mkgrmf, while effective area files 
were obtained through the fullgarf script. The spectra and 
detector responses from the ±1 orders were combined with 
the appropriate weights. Since the background is negligible, 
no correction was performed, and no systematic uncertainty 
was introduced. 

After verifying that the two observations are consistent 
with each other, in order to improve the data quality for the 
subsequent analysis we also merged the spectra (and auxil¬ 
iary files) from ObsIDs 13920/1 into a single one. It is worth 
noting that we adopted a narrower extraction strip than the 
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Figure 1. X-ray image of the central 1' X 1' region of NGC 1365 in 
the 0.3-3 keV energy range (from ObsIDs 13920/1), showing the 
extent of the diffuse starburst emission and the aperture adopted 
(green circle, 10” radius). An adaptive smoothing has been ap¬ 
plied following the same criteria of Wang et al. (2009). 

default one (20” instead of 35”): this allowed us to exclude 
part of the diffuse X-ray emission associated with the star¬ 
forming ring of ~15” diameter on kpc scale (1” ~ 90 pc. 
Fig. 0 see also Wang et al. 2009), and at the same time to 
extend the final HEG spectrum up to ~8 keV. In all the 
spectral fits (carried out with xspec v12.8, with solar abun¬ 
dances from Wilms, Allen & McGray 2000) we considered 
the 0.5-5 keV (MEG) and 1.5-8 keV (HEG) energy ranges, 
over which a total of ~1580 and 1050 counts were collected, 
respectively. We initially binned the data to 2048 channels, 
corresponding to a A A of 20 (MEG) and 10 mA (HEG), 
that is roughly the full width at half-maximum (FWHM) 
resolution. However, due to the limited number of counts, 
we further imposed a minimum of four counts per energy 
bin, and made use of the C-statistic (Gash 1979). Unless 
otherwise stated, the lines’ energies are given in the rest 
frame of NGC 1365, and all the uncertainties correspond to 
the 90% confidence level for the single parameter of interest 
(AC = 2.71). For simplicity, throughout this work we as¬ 
sumed a standard concordance cosmology with Hq = 70 
km s“^ Mpc“^, Qm = 0.27, and Qa = 0.73, although 
the luminosity distance entailed for NGC 1365 (~21 Mpc) 
slightly differs from the one based on Cepheid variables of 
D = 18.6 ± 1.9 Mpc (Madore et al. 1998). 

3 SPECTRAL ANALYSIS 

The HETG spectra of NGC 1365 are shown in Fig. 
compared to the brightest flux state ever reached by the 
source during an XMM-Newton observation performed only 
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Figure 2. Chandra HETG spectra of NGC 1365, plotted in black 
(MEG) and red (HEG), respectively. The historical extremes of 
transmission- and reflection-dominated states are also shown for 
comparison. The source has the typical Compton-thick spectral 
shape at 2-10 keV, and is ~15 times dimmer than during the 
final part of the January 2013 XMM-Newton observation (blue), 
when the neutral column density dropped to TVjj ~ 10^^ cm“^. 
Indeed, the flux lies above the minimum reached in the 2007 
XMM-Newton monitoring (green) just by a factor of ~1.8. Note, 
however, that the narrower extraction region adopted for Chan¬ 
dra allowed us to disentangle the nuclear soft X-ray component 
from most of the circumnuclear starburst emission. (All the data 
were rebinned for display purposes). 

nine months later in January 2013, when the column den¬ 
sity towards the AGN was just about Ah ~ 10^^ cm“^ 
(Braito et al. 2014; Rivers et al. 2015). Chandra clearly 
met a Compton-thick state, exhibiting a faint and very hard 
(r = —1.5) continuum above 3 keV, which closely resembles 
the lowest reflection-dominated spectrum of the first XMM- 
Newton orbit in the 2007 five-day monitoring (Fig.j^ Risaliti 
et al. 2009a). In the following analysis, we first focus on the 
overall shape of the continuum, and then take into account 
the rich emission-line spectrum in the soft X-rays. 

3.1 Baseline continuum 

As a starting point we adopted a model that simply consists 
of three different AGN-related components: the primary 
power-law continuum, attenuated by a foreground screen of 
cold gas at rest in the source frame (z = 0.005457); reflec¬ 
tion off the irradiated matter, described as a pexrav tem¬ 
plate (Magdziarz & Zdziarski 1995) with solar abundances, 
inclination of 45°, and variable strength; and a soft X-ray 
power law, which is required by the modest spectral rise be¬ 
low ~1.5 keV and is intrepreted as the faint, scattered AGN 
continuum. The latter component can be fitted in several 
alternative ways, for instance as a thermal bremsstrahlung 
of kT ~ 0.8 keV, possibly associated with starburst-driven 
shocks. However, due to its negligible impact on the final 
results and to the small (~l-2 per cent) scattering effi¬ 
ciency implied, we retained for simplicity the link of this 
soft continuum with the AGN, thus assuming a common 
photon index for the various forms of reprocessed AGN emis- 
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Figure 3. Fit residuals for the different models A-E (error bars 
are omitted for clarity). Once the AGN-related continuum compo¬ 
nents (i.e. transmitted, reflected and scattered) and the 6.4-keV 
iron K line are accounted for, a wealth of narrow emission features 
is left below ~4 keV (A). The addition of a thermal component 
from collisionally ionized gas returns an appreciable improvement 
around 1 keV (B), but any further thermal contribution is not 
strongly required (G,D). An adequate fit is obtained only after 
including individual Gaussian profiles to describe the 22 remain¬ 
ing emission lines detected with a statistical significance of at 
least 95 per cent (E). 


sion (i.e. transmitted, reflected and scattered). This model 
(hereafter tagged as model A for brevity) and all the subse¬ 
quent ones also include the Galactic column density (fixed 
to Ah = 1.34 X 10^° cm“^; Kalberla et al. 2005) and the 
conspicuous iron fluorescence line around 6.4 keV, whose 
properties remain consistent throughout the htting steps. 

Model A already provides a good representation of the 
observed continuum, while the intrinsic power law is largely 
unconstrained at this stage, also because of the decreasing 
data quality at higher energies; the photon index and ab¬ 
sorbing column were then frozen to the reasonable values 
of r = 2 and Ah = 2 x 10^'* cm“^. The overall statistics 
of Cjv = 819/489 are still rather poor without accounting 
for the emission lines (Fig. [^. To check whether some of 
these can be ascribed to optically-thin, collisionally ionized 
plasma, which is indeed a major ingredient of the XMM- 
Newton/KGS spectrum (Guainazzi et al. 2009), we added 
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Table 1. Best-fit continuum parameters for the different models (A-E) and components: primary (1), scattered (2) and reflected (3) AGN 
emission, hot (4) and warm (5) thermal emission. Nn: foreground column density. F: power-law photon index. A: power-law normalization 
at 1 keV. R: reflection strength. kT: plasma temperature. Z: elemental abundance. EM: emission measure, (f) and (t) denote fixed and 
tied parameters. Each model also includes Galactic absorption and iron Ka emission. 


Component 

Parameter 

Model A 

Model B 

Model C 

Model D 

Model E 

(1) powerlawi 






Ah (10^^ cm“^) 

200(f) 

200(f) 

200(f) 

200(f) 

140l« 

r 

2.0(f) 

2.0(f) 

2.0(f) 

2.0(f) 

2-l0l°;°® 

A (10“^ cm“^ keV“^) 

■7'7+22 
‘ ' -20 

60lm 


45+20 

^0_i8 

44t® 

(2) powerlaw 2 






r 

2.0(t) 

2.0(t) 

2.0(t) 

- 

2.10(t) 

A (10~^ cm“^ keV“^) 

l.OliO.06 

0.71±0.08 

U.DO_o 08 

- 

0 54+0-06 
0'-^^-0.05 

(3) pexrav 






r 

2.0(t) 

2.0(t) 

2.0(t) 

2.0(t) 

2.10(t) 

R 

0.50+°;®° 

'^-0.21 

0.78t°;" 

1 24+0'S^ 

1 01"’'0-10 
i-O'i-0.13 

(4) apeci 






Ah (10^^ cm“^) 

- 

- 

- 

< 0.08 

- 

kT (keV) 

- 

o 

1 + 
o o 
b b 

to 

0-96«:°^ 

1 20’’"*^'^*^ 
-‘-•^'^-0.11 

0 QQ+O-OS 

^•yy_o.o9 

Zo (Zq) 

- 

1.0(f) 

1.0(f) 

0 0^+0 04 

1.0(f) 

Znb (Zq) 

- 

1.0(f) 

1.0(f) 

0 17+0-33 

O'--*-' -0.13 

1.0(f) 

Zmg (Zq) 

- 

1.0(f) 

1.0(f) 

0 36+0-22 
0'-'^o)_o 17 

1.0(f) 

Zsi (Zq) 

- 

1.0(f) 

1.0(f) 

n 4Q+0-1® 

1.0(f) 

Zpe (Zq) 

- 

1.0(f) 

1.0(f) 

0 n^'*”0 03 

0*-Uo_0 05 

1.0(f) 

EM (10®3 cm-3) 

- 

0.25±0.06 

^•^'^-0.05 

1 Sfi+0-23 

-L-ob_o.2l 

p 

00 

1 + 

o o 

b b 

Cn 

(5) ap6C2 






Ah (10^^ cm“^) 

- 

- 

- 

< 0.35 

- 

kT (keV) 

- 

- 

n 14+0-04 

U.14_0.03 

n 15 + 0-10 

U.lO_0.05 

- 

EM (10®® cm“®) 

- 

- 

Q crc:r+0.46 
'-'•^^-0.32 

> 5.7“ 

- 

Fit statistic Cfv 

819/489 

730/487 

706/485 

671/479 

442/434 


Notes: “Due to the fact that the peak temperature of this component falls outside the fitting range, the emission measure is largely 
degenerate with the column density, so that no sensible upper limit can be obtained. 


an apec thermal component (Smith et al. 2001) with so¬ 
lar abundances. In this scenario, the ionization is driven by 
the impact of free electrons with a temperature comparable 
with the energy of the observed lines. The new fit (model B) 
brings a substantial improvement, with AC = —89 for the 
loss of two degrees of freedom only, taken by the plasma tem¬ 
perature kT ~ 0.95 keV and emission measure^ The most 
affected is the ~0.6-1.2 keV (10-20 A) range, where there are 
hints of iron L-shell emission, mainly due to 3d —>■ 2p tran¬ 
sitions from Fexvill-xxill. This hot gas can be naturally 
identified with the residual contamination from the circum- 
nuclear star formation (Fig. [^. The related continuum is 
much fainter than in the large-aperture RGS data (see also 
Fig.§, and the contribution to the most prominent fea¬ 
tures is no larger than ~30-40 per cent in the case of Lya 
lines from Ovill, Nex, and Mgxil. We therefore introduced 
a second apec component of different temperature, but this 
is remarkably cold {kT ~ 0.15 keV; Table and simply 
imparts minor changes to the noisy lower end of the MEG 
spectrum {C/u = 706/485; model G). Apparent structures 

^ The emission measure (in cm“®) is defined as EM = f rieUndF, 
where and njj are the electron and hydrogen densities. 


compatible with the He-like triplets from Mgxi, Sixill and 
Arxvil are left unexplained, but no other zone of thermal 
emission is statistically required. 

As an ultimate attempt to explore the significance of 
collisionally ionized plasma, we also allowed for variable 
abundances (tied between the two components) for all the 
key a-elements (O, Ne, Mg and Si) plus Fe, and for the 
possible column density enshrouding each star-forming re¬ 
gion. The scattered AGN continuum was neglected to avoid 
any degeneracies (model D). The same set of assumptions 
was successfully applied by Guainazzi et al. (2009) in re¬ 
producing both the RGS and the previous Chandra/ACIS 
spectra of NGG 1365. With respect to their description, in 
the HETG data the separation between the plasma tem¬ 
peratures is much more pronounced, ~0.15-1.2 keV against 
~0.3-0.7 keV, while absolute and relative abundances are 
rather similar (once the conversion between the adopted 
solar standards is applied). In any case, the global statis¬ 
tics only undergo a scarce decrease (AC = —35 for Au — 
—6), and the fit stays far from being satisfactory (C/iz = 
671/479). The full details are listed in Table Equivalent 
results were achieved by letting vary a common a-element 
metallicity over the two regions {Za ~ 0.5 and 0.03). 
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Table 2. Statistically significant emission lines in the Chandra HETG spectra, after accounting for the collisional component. The most 
plausible identifications are also listed, (f), (i), and (r) refer to forbidden, intercombination and resonance transitions for He-like species, 
(u) stands for unresolved (tr = 2 eV is then assumed). Note that the equivalent widths (EW) of broad features should not be taken at 
face value, since in most cases these are blends of two (or more) narrow lines (see the text for a thorough discussion). 



-^obs 

(eV) 

(eV) 

Flux 

(10~® cm“^) 

EW 

(eV) 

Transition 

-^lab 

(eV) 

AC/Aiz 

1. 

555/:g 

21II® 

03 + 37 

410+270 

^10-220 

Ov Kq 

0 VII Hen (f) 

554“ 

561® 

54/3 

2. 

618/12 

2(u) 

111“ 

30+39 

•^^-33 

Ov K/3 

621“ 

10/2 

3. 

655I3 

2(u) 

7-olI:? 

37±34 

Oviii Lya 

654® 

6/2 

4. 

778/12 

2(u) 

2.911® 

18±15 

0 VIII Ly/J 

Fe XVIII 3s — I- 2p 

775® 

775“ 

9/2 

5. 

8191/1 

2(u) 

0 y + 2.8 
‘ -2.1 

22±18 

0 VIII Ly7 

Fexix 3s ^ 2p 

817® 

822“ 

15/2 

6. 

874+1 

6+9 

'^-3 

4-211? 

26+23 

^0-18 

0 VIII RRC 

871® 

12/3 






Fe XVIII 3d — I- 2p 

873“ 


7. 

1024±4 

'^-3 

2.911:® 

94+22 

^^-17 

Nex Lyo 

1022® 

14/3 

8. 

1096lg 

2(u) 

u.y-o 7 

11 + 10 

Neix He/3 (r) 

1074® 

6/2 

9. 

1205±6 

8ll 

1 4+1-9 
-*-•^-0.8 

24±17 

Nex Ly/3 

1211® 

11/3 

10. 

I334I2 

2(u) 

1 4+0-9 
-*-•^-0.6 

26±15 

Mg XI Hea (f) 

1331® 

24/2 

11. 

1356±3 

2(u) 

'-'•0-0.5 

IOI 7 

Mg XI Hea (r) 

1352® 

6/2 

12. 

1476±1 

2(u) 

1 9+0.5 

l--^-0.4 

29111 

Mg XII Lya 

1473® 

42/2 

13. 

1756lg 

2(u) 

0.4±2 

ioll° 

Mg XII Ly/3 

1745® 

9/2 

14. 

18641® 

241^ 

2.8±0.7 

^o_36 

Si XIII Hea (r) 

1865® 

62/3 

15. 

2 OO 6 I 3 

2(u) 

0.5±0.3 

1911? 

Si XIV Lya 

2006® 

10/2 

16. 

2187lg 

2(u) 

Q Q + 0-5 
^•^-0.6 

40l« 

Si XIII He/3 (r) 

2183® 

9/2 

17. 

24751® 

2(u) 

'-'•^-0.5 

4llS 

S XV Hea (r) 

2461® 

11/2 

18. 

2630±20 

97+26 

-19 

^•'9-1.0 

100l“° 

S XVI Lya 

2623® 

15/3 

19. 

313811? 

10+14 

-1^-10 

1 7+0.8 
1--1-0.7 

oq+110 

oc5_42 

Ar XVII Hea (r) 

3140® 

23/3 

20. 

3310II/ 

ir^+25 

"^•'9-0.9 

1051??® 

Ar XVIII Lya 

3323® 

15/2 

21. 

3873I3® 

00+67 

^'^—22 

9 q+0.6 

^■®-0.9 

-LOU_i4q 

Caxix Hea (f) 
Caxix Hea (i) 
Caxix Hea (r) 

3861® 

3883® 

3902® 

17/3 

22. 

640411? 

24+1’^’ 

^^-12 

7 0+2-8 

' •0-2.4 

2051?? 

Fei Ka 

6403® 

38/3 


References: “Kaastra et al. (2004); *’NIST (http://physics.nist.gov/asd); '^CHIANTI (Landi et al. 2012). 


Incidentally, we note that model D actually delivers a 
fairly good fit {C/v = 513/488) of the zeroth-order spectrum 
at CCD resolution]^ although with slightly adjusted values 
(kT ~ 0.9 and 0.3 keV, A^h ~ 1.5 and < 0.2 x 10^^ cm“^; 
abundances are poorly constrained but are consistent with 
being solar). We conclude that the copious line emission 
revealed by the Chandra/'R'ETG observations is dominated 
by photoionized gas directly exposed to the radiation field 
of the AGN, which is clearly brought out for the first time. 

3.2 Emission lines 

While several star-forming regions are present within in the 
inner kpc, their weak soft X-ray emission precludes an accu¬ 
rate determination of the full range of properties (e.g. tem- 

^ A merged, low-resolution spectrum was extracted from the 
undispersed images, but is not discussed any further in this work. 


peratures and abundances) of the hot diffuse gas. We there¬ 
fore built on model B, with a single thermal component of 
solar metallicity besides the transmitted, reflected and scat¬ 
tered AGN continua. A Gaussian profile was then added for 
any residual emission line, provided that each inclusion re¬ 
sulted in a AC = —6.0 or —7.8, corresponding to the 95 
per cent confidence level for either two or three parameters 
of interest (model E). Following this criterion, which is less 
stringent than usual but is more appropriate to the signal to 
noise of the current data, 22 lines are detected (Table [^. At 
a first glance, based on the most likely identifications, all the 
main transitions of a-elements from O to Ca appear to be 
present. Most lines are narrow, and some of the broad ones 
are possible blends. Whenever the width turned out to be 
consistent with zero at the 90 per cent level, cr = 2 eV was 
assumed. This does not automatically imply that a line is 
unresolved, due to its low number of counts and significance. 
We will address this point again in the next Section. 

Interestingly, all the continuum components were left 
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free to vary and returned very sensible values (Table [^, 
leading to an excellent fit with Cjv = 442/434. The photon 
index of the intrinsic AGN power-law emission is T ~ 2.10, 
and the gas column density is A^h — 1-4 x 10^"^ cm“^; the 
absorption-corrected 2-10 keV flux is of the order of ~10“^^ 
erg s“^ cm~^, in keeping with the historical behaviour of 
NGC 1365. The strength of the reflection component is or¬ 
dinary as well, i? ~ 1 , and only ~1 per cent of the primary 
AGN continuum is scattered back into our line of sight. As 
for the collisionally ionized plasma, its characteristic tem¬ 
perature converges to kT ~ 1 keV, contributing to most of 
the Fe L-shell forest and to some extent also to the Lya lines 
from the main Q-elements (see above). With this respect, the 
spectral parameters reported in Table only refer to the 
photoionized component. Below we discuss in greater detail 
the properties of the various lines and their implications on 
the physical conditions of the emitting gas. 


4 DISCUSSION 

When seen at high resolution, the soft X-ray spectra of 
obscured active galaxies display a wealth of narrow emis¬ 
sion lines with little or no continuum, arising from low- 
temperature {T < 10® K, or a few eV) gas photoionized and 
photoexcited by the primary nuclear radiation (e.g. Sako 
et al. 2000; Kinkhabwala et al. 2002; Schurch et al. 2004). 
Chandra images have also revealed the substantial spatial 
extension of the soft X-ray emission, and its striking cor¬ 
relation with the size and morphology of the Narrow Line 
Region (NLR) as mapped by HST in the [Oiii] A5007 line 
(Ogle et al. 2000; Bianchi, Guainazzi & Ghiaberge 2006). 
This suggests that the same gas shines in both [O ill] and 
soft X-rays. In general, the latter can also be a signature of 
mechanical and shock heating in star-forming regions, yet 
the contribution from collisionally ionized plasma is usually 
minor (Brinkman et al. 2002). Although detailed spectral 
diagnostics are available for the handful of brightest sources 
only (e.g. Kallman et al. 2014), this trend seems to apply 
to the entire population of local type 1.5-2 Seyfert galax¬ 
ies, with just a few exceptions (Guainazzi & Bianchi 2007). 
NGG 1365 stands out among the peculiar objects by virtue 
of the foremost influence of thermal plasma on its RGS spec¬ 
trum. According to Guainazzi et al. (2009), evidence for 
photoionized matter solely relies on excess emission in the 
forbidden (f) and intercombination (i) components of Ovil 
Hea, in N vi Hea(f), and in the tentative Ovil and Ovill 
radiative recombination continua (RRG). 

The new Chandra/’SETSG observations take advantage 
of a much better angular resolution, and allowed us to ignore 
the bulk of the circumnuclear starburst. The resultant spec¬ 
tra cannot be explained through a multi-temperature gas 
in collisional equilibrium, whose importance is conversely 
highly diminished. The emergence of the diluted emission 
from a photoionized region closer to the hard X-ray source 
tends to reconcile the nature of NGC 1365 with the av¬ 
erage properties of the local obscured AGN. In order to 
move towards a self-consistent description of the whole se¬ 
quence of lines in Table we removed all the Gaussian 
profiles from model E (except for Fei Ka) and generated 
several sets of photoionization grids with xstar (Kallman 
& Bautista 2001). We adopted an optical to X-ray spec¬ 


tral energy distribution similar to the one used in Guainazzi 
et al. (2009), based on an Qox = 1.5 energy index[^ and 
on a 1-1000 Rydberg (13.6 eV to 13.6 keV) luminosity of 
~10"^"^ erg s“^. For each turbulence velocity broadening of 
500, 1000, 1500 and 2000 km s“^ (Ic), tke grids cover five 
orders of magnitude in both column density and ionization 
parameter]^ with 5 x lO'^® < Nn/cva~^ < 5 x 10^^ and 
0 < log (^/ergcms“^) < 5, respectively. We also tested dif¬ 
ferent gas densities, typical of BLR and NLR clouds, with 
indistinguishable results. 

The inclusion of any photoionized component was ac¬ 
cepted as significant after a statistical improvement of AC = 
— 13.3, corresponding to the 99 per cent confidence level for 
four parameters of interest, i.e. Ah, normalization and 
outflow velocity. As expected, no exhaustive solution was 
found in this effort, which should be merely taken as a 
rudimentary proof of the existence of a fairly large range 
of ionization states. This notwithstanding, a passable fit 
{C/v = 605/479) was achieved with three emission phases 
with log (^/ergcms“^) > 4.1, ~3.1±0.3, and < 0.6. The 
other gas properties are not well constrained (especially the 
column densities), so we refrain from discussing their face 
values. Anyway, there is evidence for motions along the line 
of sight: two XSTAR components have a modest outflow ve¬ 
locity of Uout = 850 I 33 Q km s while the third one prefers 
an inflow at ~3500 km s“^ to a systemic velocity (AC = 7). 
This is driven by the broad complex centred at E ~ 555 eV 
(Tableand its identification with Ovil Hea, which is not 
fully reliable (see below). In this model, once the threefold 
AGN continuum is excluded, the contribution from collision- 
ally ionized plasma to the 0.5-4 keV flux is about 30 per 
cent. Overall, the main shortcoming is that the intensity of 
many of the former Gaussian lines is slightly underpredicted, 
possibly due to higher metallicity. A 3x solar Fe abundance 
in the nuclear region of NGC 1365 is indeed suggested by the 
X-ray reflection features (Risaliti et al. 2009b; Brenneman 
et al. 2013). 

4.1 Line properties 

For the sake of a qualitative assessment, we subsequently 
considered both the MEG and HEG spectra at 1024 and 
2048 channels, preserving the constant resolution in wave¬ 
length space afforded by the gratings. This choice is not 
suitable for global statistical purposes, since no photons are 
collected in several bins, especially at very low or high en¬ 
ergies. Moreover, because of the low number of counts, the 
aspect and position of some lines is contingent on the bin¬ 
ning criterion adopted. While some degree of speculation is 
admitted in this heuristic approach, we always checked the 
raw data in wavelength units to avoid exotic interpretations, 
sticking to the most conservative explanation. None the less, 
in many cases this helped us extracting more detailed infor¬ 
mation on the properties of the emission lines, for instance 
in terms of widths, velocity shifts and relative weight of the 

® The spectral index Oox = —0.384 log(L2kev/T25ooA) 

usual indicator of the optical to X-ray luminosity ratio in AGN. 

The ionization parameter is defined as ^ = Lion/ur^, where n 
is the gas electron density and r is its distance from the source 
with 1—1000 Ry luminosity Lion. 
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Wavelength (A) 

22.543 20.664 19.074 17.712 



Wavelength (A) 

9.537 9.184 8.856 8.551 8.266 



Wavelength (A) 

4.133 3.875 3.647 3.444 3.263 3.100 



Wavelength (A) 

12.398 11.808 11.271 10.781 10.332 9.919 



Wavelength (A) 

7.293 6.888 6.525 6.199 5.904 



Wavelength (A) 

2.254 2.066 1.907 1.771 1.653 



Figure 4. Enlarged view of the MEG (a-c) and HEG (d—f) spectra in the main emission bands. The thermal and AGN continuum 
components were assumed from model E, while the lines were refitted (red curve). The data were rebinned to 1024 channels, corresponding 
to ~2 times the native FWHM resolution. For the sake of a visual inspection, a 5-channel wide triangular smoothing has been applied 
to the ungrouped spectra (turquoise curve). The vertical lines mark the expected energy of the stronger transitions, (a) Due to the poor 
signal to noise, the sampling in the oxygen band is further degraded to a factor of 8, with a smoothing window of 13 channels, (b) Both 
Nex Lyo and Ly/3 are detected with a possible velocity broadening (FWHM ~ 1300 km s“^). The thermal component is not able to fully 
account for all the putative Fe lines, (c) Despite the triplet-like appearance of the emission complex at ~1.31-1.36 keV, the lower-energy 
line is not associated with Mgxi Hea and is most likely due to Fexxi 4d 2p. (d) The broad resonance line of Sixiii can be resolved into 
a blueshifted (uout ~ 1600 km s“^) triplet, (e) The double peak across Arxviii Lyo; would imply infiows/outfiows at several thousands 
km s~^, but the red horn is not confirmed in the MEG data, (f) Besides the resolved Fe Ka feature, both Fe K/3 and Ni Ka are clearly 
present although statistically not significant due to the low effective area over this range. (See the text for more details). 
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different components in He-like triplets. To this aim, we also 
created smoothed spectral profiles by applying a triangular 
moving average to the unbinned (2048 channels) data. For 
each band of interest, the results of this procedure are shown 
in Fig.|^ where the baseline continuum is drawn from model 
E and the key features have been fitted afresh. In the fol¬ 
lowing we examine every atomic species separately, and then 
try to infer a more general picture. 

4-1.1 Oxygen 

Unfortunately, at the energies of the Ovil-vill emission the 
MEG spectrum is extremely noisy, and perhaps also affected 
by some background. In this case, we had to rebin down to 
256 channels and perform a comparably heavy smoothing 
(Fig.|^). Even so, it is clear that the broad feature at ~555 
eV in model E (which could be blindly associated with O v 
Ka from its centroid) presents several substructures, and 
that the line at E ~ 618 eV in Table (in formal agree¬ 
ment with O V K/3) is actually a blend (possibly including 
N VII Ly 7 ). We attempted at fitting the former complex with 
a self-consistent Ovil Hea triplet, fixing the energy of the 
three components and allowing for a common shift. As an¬ 
ticipated above, this selects a redshift with velocity of ~5500 
km s“^, which would be loosely supported by the identifi¬ 
cation of the line at i5 ~ 581 eV (falling just below our 
detection threshold) with Nvil Ly/3 {E ~ 593 eV at rest). 
Bearing in mind some claims in this sense (like for Mrk 110; 
Boiler, Balestra & Kollatschny 2007), such a fast inflow is 
not regarded as genuine in this context, and the ~555-eV 
blend remains ambiguous. We note that one of the spikes 
in the smoothed spectrum has the right energy for the for¬ 
bidden component of the oxygen He-like triplet, but other 
species might be involved (for example, N vi RRC is found 
at 552 eV). The statistics are seriously deficient, but the fea¬ 
ture at E ~ 581 eV can be tentatively ascribed to outflowing 
O VII. Also the O viii Lya line suggests a Vout ~ 500 ± 200 
km s“^, and the higher order Ly/3, Ly 7 and RRC transitions 
of this ion are consistent with this valuej^even if suffering 
from contamination by the Fe L-shell array. 

4-1.2 Neon 

There is very little evidence for Neix emission, once the al¬ 
most coincident 3d —>■ 2p lines from Fexix are attributed 
to collisionally ionized gas. Nex Lya and Ly/3, instead, 
are both obvious and definitely interesting. They appear to 
be resolved (a = I.OIq 5 eV for AC = 1) and minimally 
blueshifted (Uout < 400 km s“^; Table |^. The intensity ra¬ 
tio Ly(?/LyQ is quite large, ~0.6±0.3. This is a hint of pho¬ 
toexcitation, which boosts the higher order lines through ra¬ 
diative decay (e.g. Kinkhabwala et al. 2002). Also Fe L-shell 
resonances would be enhanced (Sako et al. 2000), and indeed 
the thermal component fails to reproduce their emission at 
~ 1.06-1.18 keV (Fig. Notably, the strongest 3d —> 2p 
Fe lines, including Fexx (0.956 keV) and Fexxi (1.009 keV), 
accept some blueshift at Vout ~ 900 km s“^. This is proba¬ 
bly too high (e.g. Veilleux, Cecil & Bland-Hawthorn 2005) 

® In the Ovill energy range, the MEG absolute wavelength ac¬ 
curacy is equivalent to >^200 km s“^. 



Rest-frame energy (keV) 

Figure 5. Sixill Hea triplet in the HEG spectrum at 2048 chan¬ 
nels (AA ~ 10 mA, consistent with the FWHM resolution). The 
blue and green dashed lines indicate the energy of the forbid¬ 
den, intercombination, and resonance transitions for Uout = 0 
and 1600 km s“^, respectively. Both velocity components seem 
to be present within the energy/wavelength calibration accuracy, 
the former possibly associated with the starburst and the latter 
to the AGN. The continuum is drawn in red. 

for any starburst superwind in NGC 1365, favouring an ori¬ 
gin from matter photoionized and photoexcited by the AGN 
rather than shock-heated much farther out. 

4 . 1.3 Magnesium 

Three lines are seen in the Mgxi Hea band (Fig.|^), but the 
low-energy one is actually due to Fe xxi 4d —>■ 2p, again more 
prominent than in a pure thermal scenario. The remaining 
pair are presumably the forbidden and recombination com¬ 
ponents of the triplet. Together with the Mgxil Lya line, 
they share a — 1.6 ± 0.6 eV (AC = 1) and Uout ~ 700 ± 100 
km s“^ (Table B. The only other feature from these ions is 
the possible Mgxil RRC, while the puzzling excess around 
the energy of the Ly/3 is formally incompatible with the shift 
of the parent series, and is eventually unknown (Fig. [^1). 

4 . 1.4 Silicon 

A broad line is detected in model E at the energy of Sixill 
Hea, centred around the resonance transition. Its width in 
the HEG spectrum at 1024 channels is cr ~ 18 eV, and a 
narrow core is noticeable (Fig. [^1). While the FWHM of 
~7000 km s“^ would dictate a BLR origin, we tried to dis¬ 
entangle the triplet components. The best solution is found 
for unresolved profiles at Vout ~ 1600±100 km s“^, with the 
forbidden and resonance lines shifted to ~1.850 and 1.875 
keV (Fig.^. This outflow velocity is not confirmed by Sixiv 
Lya, which lies at the expected rest-frame energy. Moreover, 
after the joint ~10-eV shift, the intercombination line be¬ 
comes the stronger one at ~ 1.865 keV, which is inconsistent 
with the weights of Mg xi and practically unphysical. On the 
other hand, its intensity is naturally explained by a contribu¬ 
tion from resonant emission at rest, supported by the faint 
forbidden companion and associated with either the AGN 


© 2015 RAS, MNRAS 000.[lHTT] 













































Chandra/HETG view of NGC1365 9 


Table 3. Summary of the best constrained lines. No errors are 
given for tied (or frozen) parameters. Uncertainties are purely 
statistical at the 68 per cent confidence level (AC = 1). 


Line 

-^lab 

<7 

FWHM 

■^OUt 

(keV) 

(eV) 

(km s“^) 

(km s“^) 

Ne X Lya 

1.022 

1 9+0-7 
-*-•^-0.5 

1280t«° 

< 400 

Nex Ly/3 

1.211 

2.1 

1280 

< 400 

Mg XI Hea (f) 

1.331 

1.6±0.6 

o 3 n-|-c5UU 

o4U_32o 

710±110 

Mg XI Heo (r) 

1.352 

1.6 

840 

710 

Mg XII Lya 

1.473 

1.8 

840 

710 

Si XIII Hea (f) 

1.839 

< 1.6 

< 600 

15901®° 

Si XIII Hea (r) 

1.865 

< 1.6 

< 600 

1590 

Si XIV Lya 

2.006 

< 1.7 

< 600 

< 300 


or the starburst. Two velocity components are tentatively 
distinguished then. Only the properties of the blueshifted 
triplet are listed in Table Finally, there is an unidenti¬ 
fied secondary peak at if ~ 2.039 keV. The closer transition 
is Alxill Ly/3 {E ~ 2.048 keV), but any Lya counterpart 
would be much weaker. If due to H-like silicon, the sizable 
shift calls for Uout ~ 4700 km s“^. 

4-1.5 Argon 

While S xv-xvi is recognized but very faint, some basic con¬ 
straints can be put on Arxvil-xvill. The broad (cr ~ 17 eV) 
Hea line is found again at the resonance energy. The same 
considerations made above for the Si triplet could apply, but 
here any fine structure cannot be appreciated since the sep¬ 
aration of the three components is commensurate with the 
FWHM. The same holds for He/3. Another double-peaked 
feature is visible across the Arxvill Lya transition (Fig.[^). 
The red horn is dubious, since the HEG and MEG data are 
in conflict (the MEG resolution, however, is lower by a factor 
of 2). The blue one would correspond to Pout ~ 2700 km s“^. 
In any case, we can rule out a contribution of the Arxvill 
Ly/3 line (3.936 keV at rest) to the blend near ~3.90 keV, 
which is the alleged combination of forbidden and resonant 
Caxix Hea with modest blueshift. 


4.1.6 Iron 

Above 4 keV, only the fluorescence Ka line from neutral iron 
at if = 6.404^0 015 is statistically significant, although the 
K/3 is clear as well, and there is a trace of Ni Ka (Fig. |^). 
All of these are signatures of reflection of the AGN radia- 
dion in cold, optically thick matter, arguably located at the 
torus scale in the original unification scheme. We can as¬ 
sess the distance of this gas from the width of the Fe Ka 
line, which is nicely resolved thanks to its ~35-40 counts 
(cr = 27^12 eV; Fig. §. We assume that the broadening is 
due to the virial motion in the emitting region, and that the 
velocity dispersion is related to the FWHM through a factor 
/ = %/3/2, which corrects for the (unknown) geometry and 
distribution of the gas (see also Netzer & Marziani 2010). 
The predicted FWHM is ~760(M8/rpc) km s“^, where Mg 
is the mass of the central black hole in units of 10® Mq and 
Tpc is the radial distance in pc. This expression agrees within 
a factor of ~2 with the case of pure Keplerian motion with 



FWHM (km s-^) 


Figure 6. Confidence contours at the nominal 68, 90 and 99 per 
cent level for the Fe Ka intensity versus velocity width, obtained 
from the unbinned HEG spectrum with all the AGN continuum 
parameters free to vary. The vertical dashed line indicates the 
FWHM of the broad H/3 component in the optical, suggesting a 
common origin from BLR clouds. 

reasonable orbital inclination {i > 20°). Far larger uncer¬ 
tainties are associated with the black hole mass. No direct 
measurement (e.g. from reverberation mapping) is available 
for NGG 1365, so that the full range of log (Mbh/Mq) ~ 6- 
8 is covered in the literature. The higher values are based on 
the correlation of Mbh with the bulge luminosity, and are 
likely biased by the starburst contribution to the latter. The 
most sensible guess is then of a few x 10® Mq (Risaliti et al. 
2009b; Davis et al. 2014). With a conservative Mg = 0.1, and 
the observed Ka FWHM of ~30001J^35 q km s“^, we obtain 
that r ~ 0.002-0.02 pc, which is within the BLR. Indeed, 
the width of the Ka line is consistent (at 90 per cent errors) 
with that of the broad nuclear component of the H/3 in the 
optical, whose FWHM is 1895 km s“^ (Schulz et al. 1999). 
From the Ka/H/3 width ratio we can deduce that the relative 
size of the respective emission regions is roughly the same. 
This is not surprising in a source like NGG 1365, where BLR 
clouds are responsible for variable X-ray obscuration up to 
the Fe-K band and beyond. Regardless of the precise black 
hole mass, in terms of the gravitational radius the distance 
of the Fe Ka reflector perfectly agrees with the universal lo¬ 
cation of a few X10^ Tg determined by Shu, Yaqoob & Wang 
( 2011 ). 


4.2 Covering factor 

As established from the FWHM of Fe Ka (3000 km s“^) 
and, provisionally, from those of Nex (~1300 km s“^) and 
Mg XII Lya (~800 km s“^), the gas probed by these Chan¬ 
dra/}AETG spectra is likely located across the boundary be¬ 
tween the classical BLR and NLR (see also Grenshaw & 
Kraemer 2007). It is therefore tempting to try to evalu¬ 
ate its covering fraction (/cov) and to explore its connection 
with the clumps responsible for the X-ray occultations of the 
source. To do so, we applied an xstar grid to the 1.3-1.5 
keV band, characterized by the intense lines from Mgxi-xil. 
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For given column density and amount of ionizing radiation, 
the measured Mg luminosity of Lms ~ 6 x 10®® erg s“® de¬ 
pends on the fraction of solid angle occupied by the gas. 
The predicted luminosity for a full covering of 47r sr (i.e., 
/cov = 1) is degenerate with the column density, but be¬ 
comes slowly responsive to its exact value at A^h > 2 x 10 ®® 
cm“®, which is typical for BLR clouds in NGC 1365 (Risaliti 
et al. 2009b). In that range, the entailed /cov lies between 
~2-7 X 10“®. This is not far from the NLR covering fac¬ 
tor delivered by photoionization models, which, however, is 
thought to be substantially underestimated (Netzer & Laor 
1993). Indeed, it should be quite similar to the BLR one 
(~0.2-0.4; e.g. Korista, Ferland & Baldwin 1997), even if 
showing a stronger decreasing trend with luminosity (Stern 
& Laor 2012). A column of ~10®® cm“®, still plausible for 
the NLR, would be needed to match the expectations. 

There are several caveats about this simple calcula¬ 
tion, though. For instance, the nearly equivalent intensity 
of the Hea (f) and Lya lines requires a fine tuning of 
log (^/ergcms”®) ~ 3.5 in a single-zone model. If a wider 
range of ionization states is involved, /cov can be consid¬ 
erably higher. Allowing for two components with Ah = 

5 X 10®® cm~® and log (^/ergcms”®) = 2.5 and 4.5, respec¬ 
tively, returns a better fit and a covering fraction of ~0.15. 
Many other configurations are viable, and also a mild at¬ 
tenuation of Livig by an external absorption layer, maybe at 
host-galaxy scales, cannot be ruled out. Alternatively, given 
the narrower profile of the Mg lines compared to Fe Ka, we 
can speculate that the innermost clouds significantly reduce 
the ionizing flux seen by the emitting gas at larger distance, 
thus producing an apparently low /cov for the NLR. For 
log (^/ergcms”®) ~ 3.5 and n ~ 10® cm“®. Mg emission 
would occur at r ~ 0.2 pc, at least an order of magnitude 
higher than what inferred for Fe Ka (remarkably propor¬ 
tional to the square of the FWHM ratio). This is the most 
intriguing explanation, since it is consistent with the under¬ 
luminosity in the soft X-ray oxygen lines noted by Guainazzi 
et al. (2009), despite the fact that the AGN in NGC 1365 is 
not intrinsically weak. Unfortunately, the current data qual¬ 
ity does not warrant any conclusive answer. 

4.3 Temperature and density 

The presence of possible RRC from species like O viii or 
Mg XII implies that the emitting gas has a low tempera¬ 
ture and that the ionization state is regulated by radiation 
rather than by collisional equilibrium. Radiative recombina¬ 
tion, in fact, is the emission of a photon following the cap¬ 
ture of a free electron, mainly to the ground state of an ion. 
The shape of the consequent feature is a powerful diagnos¬ 
tic of the gas temperature. In a cold plasma, only the tail of 
high-energy electrons overcome the recombination thresh¬ 
old, and the RRC have a narrow, line-like profile, whose 
width is AE ~ kT (Hatchett, Buff & McCray 1976; Liedahl 

6 Paerels 1996). Both detections in proximity of Ovill and 
Mg XII RRC rely on ~6-7 line counts only in the MEG and 
HEG data, respectively. Hence we cannot disentangle the 
putative Ovill RRC from the nearby Fexvill L-shell tran¬ 
sition, in particular if the former is moderately blueshifted. 
No blend with other strong lines is instead foreseen at the 
energy of Mgxil RRC, yet we remind that this feature is 
not statistically significant, so we were only able to put an 


upper limit of 3 eV (or 3.5 x 10^ K; AC = 1) with a redge 
model. 

More information can be potentially derived by the 
He-like triplets, consisting of the resonance (ls2p ®Pi —>■ 
Is® ®So), intercombination (ls2p ®P 2 ,i —t Is® ®So), and for¬ 
bidden (ls2s ®Si —>■ Is® ®So) transitions. These provide an 
effective measure of the physical conditions of the gas (Por- 
quet & Dubau 2000, and references therein), since the in¬ 
tensity ratios 77. = f/i and Q = (f-|-i)/r are very sensitive 
to the electron density and temperature, respectively. Given 
that O VII is too noisy, Ne ix is formally undetected, and 
Sxvil, Arxvil and Caxix are too faint and/or unresolved, 
the only accessible species are Mg xi and Si xill. Even if the 
triplet decompositions proposed above (which impose out¬ 
flow velocities of ~700 km s“® for Mg and ~1600 km s“® 
for Si) are valid, only upper limits can be placed. The den¬ 
sity can be anything below the critical value (n ~ 10®®-10®'’ 
cm“®) beyond which the forbidden components are colli- 
sionally suppressed, while temperature is constrained to be 
T < 3 X 10® K. This is barely meaningful in the light of 
the appearance of narrow RRC features, but it is enough to 
exclude collisional ionization in a hot {kT ~ 1 keV) plasma 
as the dominant mechanism in these emission-line HETG 
spectra of NGC 1365. 


5 SUMMARY AND CONCLUSIONS 

In this paper we have presented the X-ray spectral analysis 
of the first Chandra/}TETG observations of the prototypi¬ 
cal changing-look Seyfert galaxy NGC 1365, which lingered 
in a Compton-thick state for the entire span of four days. 
The soft X-ray spectrum was thus dominated by a wealth 
of recombination lines belonging to He- and H-like ions of 
the most common light elements from oxygen to calcium, 
plus iron L-shell transitions. Thanks to the unrivaled spa¬ 
tial resolution offered by Chandra, we were able to zoom in 
on the close surroundings of the AGN, and unveil a pho- 
toionized gas component that was largely diluted in previ¬ 
ous observations. Indeed, the spectra obtained with XMM- 
Newton/KGS, whose large aperture also encompasses all the 
diffuse emission from the inner 5 kpc, were heavily contam¬ 
inated by collisionally ionized plasma, shock-heated by the 
herce circumnuclear star-formation activity. The residual 
thermal contribution to the overall emission-line intensity 
now amounts to ~30 per cent only. 

The emergence of a photoionized spectrum with prop¬ 
erties analogous to those reported in most of the nearby 
obscured Seyfert galaxies partly mitigates the apparent 
anomaly of NGG 1365. In spite of the modest statistical 
quality of the data, whereby only a handful of lines have 
> 10 counts, the picture in favour of cold {kT E) gas ex¬ 
posed to the AGN radiation is corroborated by the tentative 
presence of narrow RRC features (Ovill and Mgxil), and 
by the qualitative inspection of the He-like triplets (mainly 
from Mg and Si). Photoexcitation might be important as 
well, as indicated by the strength of some higher order tran¬ 
sitions (e.g. Nex Ly/3) and of the Fe L-shell forest around 
~1 keV. Possible outflow velocities in the range ~0-1600 
km s“® are revealed, and a few lines (Nex, Mgxi-xil) show 
some evidence of broadening (~1000 km s“®), signifying 
that the location of the gas is across the virtual BLR/NLR 
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boundary. The Ka fluorescence feature from neutral iron at 
E = 6.404^0 015 keV is resolved to a width of u = 271^2 
eV, corresponding to a FWHM of ~3000ti35o km s“^. This 
is somewhat larger but fully consistent with the optical H/3 
line, supporting a comparable size of the emitting regions. 
With a preferential distance of few x 10'* rg from the black 
hole, it is then very likely that Fe Ka arises from the same 
Compton-thick clouds that induce the recurrent column den¬ 
sity jumps and the extreme variability seen in this source. 

Besides the main findings, there are also some unusual 
hints, such as of broad (BLR-like) profiles in resonance lines, 
of double peaks, and of radial velocities of several thousands 
km s“^, even in inflow. If confirmed, these would be abso¬ 
lutely exceptional traits, but caution is mandatory given the 
poor statistics. As opposed to standard type-2 Seyferts like 
NGC 1068, the erratic behaviour of NGC 1365 below 10 keV 
makes it hard to amass strictly consistent information from 
different epochs, so naively it could seem unfeasible to over¬ 
come this limitation with just a longer exposure. However, 
even in a Compton-thin case, the column density is normally 
of the order of A^h ~ 10^® cm“^, and the direct AGN contin¬ 
uum is almost completely absorbed in the soft X-rays. The 
ample record of past observations suggests that a nearly un¬ 
obscured state such as the one sampled in January 2013 is a 
very rare occurrence for NGC 1365. The photoelectric cutoff 
is thus expected to fall beyond the main complexes from O, 
Ne, Mg, and maybe also Si, leaving their aspect virtually 
constant with time and allowing the build-up of a nuclear 
emission-line template. On the other hand, the appearance 
of the transmitted AGN continuum at higher energies would 
enable the study of the high-ionization absorption features 
originating from the accretion-disc/BLR wind, from Fe xxv- 
XXVI down to Gaxix-xx and even Arxvii-xviii, depending 
on the exact value of Ah. For these reasons, any prospec¬ 
tive observation of NGG 1365 at high spatial and spectral 
resolution would provide further invaluable insights into the 
properties of the sub-pc scale environment in active galaxies. 
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